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Abstract — This paper presents an approach based on Genetic Algorithm to soh^VtheNconomic load 
dispatch problem with losses for three thermal plant system and combined cyc«^^eneration plant. 
Genetic algorithms are adaptive search methods that simulate some of the naJsraLprocesses: selection, 
inheritance, random mutation and population dynamics. This approach is usptH^cest for an example of 
three thermal plant system and one combined cycle co-generation plant in*tj/^ways of dispatching load 
between three thermal plants and among two thermal plus one combinsd/^d»co-generation plant and the 
results are compared. In this paper, the results are obtained throughtj^glnetic algorithms developed in 
c Language. ?S> 

Index Terms — About four key words or phrases in order of imj^^knce. separated by commas. 



Economic load dispatch (ELD) is a sub-problem of A^Optimal power flow (OPF) having the objective of 
fuel-cost minimization. The classical solutions for problems have used equal incremental criterion for 
the loss-less system and the use of penalty factorWor considering the system losses. However, all these 
methods are based on the assumption of cctffir^ity and differentiability of cost functions. Hence, the cost 
functions have been approximated in the^lffe^entiable form, mostly in the quadratic form. Further, these 
methods also suffer on two main coajfcsjOne is their inability to provide global optimal solution and 
getting stuck at local optima. The sScKd^ problem is handling the integer or discrete variables. 

Genetic algorithms (GAs) have^^n proved to be effective and quite robust in solving the optimization 
problems. GAs can prowl^nifir global solutions and also can handle effectively the discrete control 
variables. Discontinuity afjcNaoli-differentiability of cost characteristics can be effectively handled by GAs. 



ation plants (CCCP) have the following advantages 



• the thermal plants, 



Combined cycle 
namely, 

gjfMl thermal efficiency. 

air pollution by nox, dust etc. 
ndent operation of gas turbine for peak loads, 
start-up and less capital cost per kw and 
water requirement per unit of electrical output. 

This paper proposes the application of GAs to solve the economic load dispatch for two types, that is, 

1. The thermal plant systems and the results are compared with conventional method, and 

2. Two thermal plant systems and third plant as a combined cycle co-generation plant and the results 
are included. 
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II. Theoretical Analysis on Genetic Algorithm and Economic Load Dispatch 
A. Genetic Algorithm 

l) Introduction 

Genetic algorithms are search algorithms based on the mechanics of natural selection and natural genetics. 
They combine survival of the fittest among string structures with a structured yet randomized information 
exchange to form a search algorithm with some of the innovative flair of human search [i]. 
generation, 



Fx 



Figure i The Single Peak function 



A new set of artificial creatures (strings) 
occasional new part is tried for good measi 




The mechanics of simple genei 
copying strings and swappini 
attractions of the genetic 



;) ispf« 



t easy for Calculus - Based Methods 

! bits and pieces of the fittest of the old; an 



imple Genetic Algorithm 



;orithm are surprisingly simple, involving nothing more complex than 
strings. Simplicity of operation and power of effect are two of the main 
im approach. 



A simple genetic ajep/ahn that yields good results in many practical problems is composed < 
operators: ^^V^ 



<3 



ction 
Nfesover 

utation 



A) Reproduction 



Reproduction is a process in which individual strings are copied according to their objective function 
values, f (biologists call this function the fittest function). This operator is artificial version of natural 
selection, a Darwinian survival of fittest among string creatures. In natural populations fitness is 
determined by a creature's ability to survive predators, pestilence and the other obstacles to adulthood and 
subsequent reproduction. In our artificial setting the objective function is the final arbiter of the string 
creature's life or death. 
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Kafs. The easiest is 
<tte wheel slot sized 
Jack box problem has 



Table i 

The reproduction operator may be implemented in an algorithmic form in a numbei 
to create a biased roulette wheel where each current string in the population has 

in proportion to its fitness. Suppose the sample population of four strings in ti 

objective or fitness function values f as shown in Tablei ^fa* 

Summing the fitness over all four strings, one obtains a total of ii7Q*jCj>ercentage of population total 
fitness is also shown in the table. To reproduce, one simply spin tft^^jrlted roulette wheel thus defined 
four times. For the example problem, string number i has a fitrfes^alue of 169, which represents 14.4 
percent of the total fitness. As a result, string 1 is given 14.4 perfon* of the biased roulette wheel, and each 
spin turns up string 1 with probability 0.144. 1 



After reproduction, simple crossover may pro^SpVinto two steps. First, members of the newly of 
reproduced strings in the matting pool are JBattea at random. Second, each pair of strings undergoes 
crossing over as follows: an integer positio^klflong the string is selected uniformly at random between 1 
and the string length less 1 (1 to l-i). Swaj^Ag a U characters between positions k+i and 1 inclusively creates 
two new strings. For example conside^rVujigs Ai and A2 from the initial population: 



Ai=ono/i 

A2=110l/0 

Suppose in choosinj 
symbol '/')• The resul 
the ne^ 



a rtnasi 



<2T 



Ai'=onoo 
A2W1011 



im number between 1 and 4, we obtain a k=4 (as indicated by a separator 
Tossover yields two new strings where the prime (') means the strings are part of 



^*BiAj^ population of n strings over some appropriate alphabet, coded so that each is a complete idea or 
^scnption for performing a particular task. Sub strings within each string (idea) contain various notions 
jf wMat is important or relevant to the task. Genetic algorithms ruthlessly exploit this wealth of information 



(1) Reproducing high-quality notions according to their performance and 

(2) Crossing these notions with many other high-performance notions from other strings. 

Thus, the action of crossover with previous reproductions speculates on new ideas constructed from the 
high-performance building blocks (notions) of past trials. 
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C) Mutation 



Mutation plays a decidedly secondary role in the operation of genetic algorithms. Mutation is needed 
because, even though reproduction and crossover effectively search and recombine extant notions, 
occasionally they may become overzealous and lose some potentially useful genetic material (i's or o's at 
particular locations). 

One may notice that the fitness or objective function values are the same as the black box values (compare 
Tables i and 2). There is no coincidence, and the black box optimization problem was well repres^nj/oNby 
the particular function f(x), and coding one is now using. A generation of the genetic algorithm begf^^ith 
reproduction. One select the mating pool of the next generation by spinning the weighted roi^t^e«wheel 
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Actual simulation of this process using coin tosses kAresulted in string 1 and string 4 receiving one copy in 
the mating pool, string 2 receiving two copies, anchoring 3 receiving no copies, as shown[9] in the center of 
Table 2. 



omic Load Dispatch 
) Introduction 

An early method of attemptiwAta/rninimize the cost of delivered power called for supplying power from 
only the most efficient pl«tJO'gh ,: loads. As load is increased, the most efficient plant would supply power 
until the point of maxijrjAmyfefftciencv of that plant was reached. Even with transmission losses neglected 
this methods fail to mfoArize the cost. 



following trends in the growth of power systems it has become progressively 
icreasing attention to economic operation of power systems. 

1. ^*kn%ny cases, economic factors and the availability of primary essential such as coal, water etc., 
^^dfttate that new generating plants be located at greater distances from the load centre. 

The installation of larger blocks of power has resulted in the necessity of transmitting power out of 
a given area until the load in that area is equal to new block of installed capacity 

3. Power systems are interconnecting for purposes of economy interchange and reduction of reserve 
capacity. 

4. In a number of areas of the country the cost of fuel is rapidly increasing. 



The main factor controlling the most desirable load allocation between various generating units is the total 
running cost. The operating cost of a thermal plant is mainly the cost of the fuel. Fuel supplies for thermal 
plants can be coal/natural gas, oil or nuclear fuel. The other costs such as cost of labour, supplies, 
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maintenance etc., being difficult to be determined and approximate, are assumed to vary as affixed 
percentage of the fuel cost. 

2) Economic Dispatch Neglecting Losses 

When transmission distances are very small and load density is very high, transmission losses may be 
neglected and the optimal dispatch of generation is achieved with all plants operating at equal incremental 
production cost. .The common for including the effect of transmission losses is to express the total 
transmission losses as quadratic function of the generator power outputs . The simplest quadratic fornyTV 

V 

n n + 

P1=Z Z PiBijPj (1) 

c \ 

Here the coefficients Bij are called loss coefficients or B-coefficients. B-coefficients\f\a turned constant 
and reasonable accuracy can be expected provided the actual operating condition^J^ close to the base 
values where the B-constants are computed. 

For a system with 'n' generating units let 

Ft=fi+f2+ +fn vN>^^ 

Where Ft is the cost function giving all the fuel for the entire syStenwtnd is the sum of the fuel cost of the 
individual units fi, f2 , f n. 

The total megawatt power input to the network from all 
PGi +PG 2+PG 3 PGn. 

Where PGi ,PG2 , PG n are the individual oj^tput^of the units injected into the network. The total Ft of 

the system is a function of all the power plaifljnukputs. The constraining equation on the minimum value of 
Ft is given by the power balance equatioa**^^ 

Pi+Pd- Z PGi=o ....*S(X 

Where Pd is the total pfrwLryceived by the loads and Pi is the transmission loss of the system. The 
objective is to obtain a frnWmum Ft for affixed system load Pd subjected to the above power balance 
constraint. ^Ak- 

The procedure foV^^ljnng such minimization problems is called the method of Lagrange's multipliers. 

The new casVIHuction F is formed by combining the fuel cost and the equality constraint in the following 
manner^ 




F^^t A ( Pi +Pd- Z PGi) (.3) 



The augmented cost function F is also called the Lagrangian and the parameter, which is called the 
Lagrange multiplier is the effective incremental fuel cost of the system when transmission line losses are 
taken into account. 
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The original problem of minimizing Ft subjected to the constraint given by Eq 3.2.4.2 is transformed into an 
unconstrained problem given by F, where it is required to minimize F with respect to and the generator 
outputs. 

n 

F/PGi= Ft/PGi + A ( Pl+Pd - Z PGi ) / PGi =0 (4) 

i=o 

Since, Pd is fixed and the fuel cost of any one unit varies only if the power output of that unit is variad^>o 
the above equation yields: 

F/PGi=Ft/PGi + A (Pl/PGi-i) =0 (5) &\ * 

for each of the generating unit outputs PGi,PG2, PGn. ^ ^V^ 1 ^ 

Because Ft depends on only PGi, the partial derivatives of Ft can be replaced by tha^LHierivative and the 
Eq 3-2-4-5 gives /^S^ 



X = (dFt/dPGi) (6) 

i-(Pl/PGi) 
for every value of i. 



This equation can also be written in the form ^TJ 



1=L i (dFt /dPGi) (.7) 

where Li is called penalty factor of plant i, and is giv 



(.{ 

i-(Pl/PGi) 



L^rrmrr 



The result for Eq (7) means that thMlmWura fuel cost is obtained when the incremental fuel cost of each 
unit multiplied by its penalty fatoor lWne same for all the generating units in the system. The products Li 
(dFt/dPGi) are each equal to c^le^i the system, which are approximately the cost in Rs/h to increase the 
total delivered load by 01^ 



Eq (7) governs the romAjjJtion of transmission loss into the problem of economic loading of units in 
plants, which are geggtj/^ically dispersed throughout the system. 



<9 



ine^^^e 

J4T 



3) Combined Cycle Co-Generation Plant 

co generation plants have the following advantages over the thermal plants, namely, 



;her overall thermal efficiency 
minimum air pollution by nox, dust etc. 

3. independent operation of gas turbine for peak loads. 

4. quick start-up and less capital cost per KW and 

5. less water requirement per unit of electrical output. 



The fuel consumption and the cost characteristics of such plants, in general, are not differentiable. The 
application of the genetic algorithms is the only viable solution for power system with combined cycle 
cogeneration plants and it is not possible to solve such a ELD problem by conventional technique. Till date 
only limited work has been reported in the area of CCCP economic load dispatch. 
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4) Two Thermal Plants and One CCCP System 

In the three thermal plant systems, the third thermal plant is to be replaced by combined cycle 
cogeneration plant CCCP( two 75MW gas turbines and one 50 MW steam turbines). The fuel cost 
characteristics of this plant is shown in figure below[6] : 




POWER (MW) V^^fc' 

)f co- ^^ jaym plant 
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Table.3 
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Using the given set of points, a polynomial of specific order can be fitted using polynomial regression 
technique. We assume that the quadratic is to be fitted here for the cost curve. Assume that n pairs of 
coordinates (xi,yi) are given which are to be approximated by a quadratic. Let the quadratic be represented 
as 



Y = AXz+BX+C 



When X=Xi The left hand side of the above equation is represented by Yci 

Yci=AXi2+BXi+C ^^V^ 
The sum of squares of the deviations is given by 

S=Z (Yi-Yci)2 = Z(Yi-AXi2-BXi-C)2 % 

Differentiating S w.r.t A,B,C respectively and setting each of the co-efficient equal to ^gy, we get 

nC+BZXi+AZXi2 = Z Yi Qr 
CZXi+BZXi2+AZXi 3 = ZXi Yi ~ 
CZXi2+BZXi 3 +AZXi 4 = ZXi 2 Yi \Jv^ 



These are three linear equations in three unknowns. These arejca^d normal equations for quadratic 
regression. These may be solved by Gauss - Jordan procedure. SSi 



Finally we obtain the curve as 

F = A PGi2 + B PGi 



ircmOT*he ger 



For the problem, the solution was obtained thrcmah^he genetic algorithms developed in c language. It is 
observed that this method is accurate and may replace effectively the conventional practices presently 
performed in different central load dispatch/PVpfVrs. 

III. Cla^^^^lonomic Load Dispatch Problem 

The object of ELD problem is toimuiimize the total fuel cost at thermal plants 

OBJ=Z Fi(Pi) C/> 



^C»Pi-Pl-Pd=o 

Therneq 



Thernequalities of real power limits on the generator outputs are 
Pimin < Pi < Pimax 



Where Fi(Pi) is the individual generation production in terms of its real power generation Pi, Pi is the 
output generation for unit i, n the number of generators in the system, Pd the total current system load 
demand, and PI the total system transmission losses. 
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The thermal plant can be expressed as input-output models (cost function), where the input is the fuel cost 
and the output the power output of each unit. In practice, the cost function could be represented by a 
quadratic function. 

Fi(Pi)=Ai.Pi 2 +Bi.Pi+Ci 

The incremental cost curve data are obtained by taking the derivative of the unit input-output equation 
resulting in the following equation for each generator: 



dFi(Pi)/dPi=2.Ai.Pi+Bi 



transmission losses are a function of the unit generations and are based on the system topoloj ine the 

ELD equation for a specified system requires an iterative approach since all unit generati*n\^j|>cation are 
embedded in the equation for each unit . 

Pl=PiBijPj /»^^ 

Where Bij are coefficients, constant for certain conditions. *fV 
Application of GA to economic load dispatch problem. * 

Value=bito*20+biti*2i+ +biti*2i+ +bitchrom-length*2chrom-l%^^ 

If the optimized parameter belongs to (Pimax,Pimin),decodin^^lue of the parameter is computed by 
equation(4.i.i). \ 4 

value*(Pimax- Pimin) » 
Pi = Pimin + 



2chrom-length -l 

Objective Functio^SyiTl Fitness Function Formulation 



In the ELD problem, the goal is to r*i^m^r the objective function 



with the constraint 



off3fogIity 



,6 



to constrained optimization problem and thus forming fitness function 
Fct=Ft+PF (EPi-PL-Pd) 

where pf is penalty factor .the penalty function is placed into the objective function in such away that it 
penalizes any violation of the constraint and forces that unconstrained optima toward the feasible region . 
in the ELD problem the goal is to minimize the objective function Fct, while the objective when using GAs 
is to maximize a fitness function Ftt in the given form[io]. 

Ftt=EXP(-(Ki*Fct)k2) 
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Ki and I<2 are constants and the value is problem dependent .considering the evolutionary process of the 
GAs ,the solution is improved through the generation and also to decrease the penalty function over 
successive iterations can be adapted with the penalty function varying directly with the number of 
generations. 



IV. Software and Results 



The cost functions of the three thermal plants considered in this paper are obtained from Sheble and 
Britting [n]and they are as follows. 

Fi=(o.ooi56)Pi A 2+(7.92)Pi+56i Rs/h 
F2=(o.ooi94)P2 A 2+(7.85)P2+3io Rs/h 
F3=(o.oo482)P3 A 2+(7.97)P3+78 Rs/h 
The loss co-efficients of the considered system are 
0.000075 0.000005 0.000008 
0.000005 0.000015 0.000010 
0.000008 0.000010 0.000045 
The Operating ranges for three plants are : 

100MW < Pi< 600MW 

100MW < P2 < 400MW 

50MW < P3 < 200MW 

Population size :io 
Chromosome length 36 
Maximum number of generations : 1000 
Crossover probability 10.50 
Mutation probability : 0.01 
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3i$;: :;■;$«: i;ss:: 



:.-:=< 



Table.4 

Here the range of limits for the load can be taken as 250MW - 1200 MW. 

Genetic Algorithm claims to provide near optimal or optimal solution for coi 
problems. Therefore the effectiveness of the genetic algorithm solutions should 
experimental results. For economic load dispatch problem, the results obtained 
developed in C language was tested for three thermal plant systems and exten| 
cycle co-generation plant in three thermal plant systems. The execution time t, 
120 MHz processor. 



m«ti#nally intensive 
"^ys be evaluated by 
genetic algorithms 
one plant as combined 
two seconds in Pentium 



The performance of GA approach is compared with the classical I^cVnayer method as given in table 5.1. 
For the load of 585.33 MW, ZPi-PD-PL = 592.297 - 585.33 - fjfifra = 0.005 MW in Classical Kirchmayer 
method, and it is 592.850 - 585.330 - 7.25 = o in GA appr^crr^For the load of 869 MW ZPi-PD-PL = 
884.390 - 869.000 - 15.420 = -0.030 in conventional methocfe^d it is 884.811 - 869.000 - 15.812 = -0.001 MW 
in GA approach. The constraint of equality ZPi-PD-HAeiydoes n °t hold good for classical Kirchmayer 
method. But this constraint is almost equal to or neA^Eo zero in GA approach. It at all it is having some 
value, it is multiplied by the penalty factor (PF) ^Nk> is placed into the objective function in such a way 
that it penalizes any violation of the constraints ancrforces that unconstrained optima towards the feasible 
region. The solution of GA is improved iJtoHgh the generations, and also by decreasing the penalty 
function. 

In conventional method, the value»^or^ > Gi,PG2,PG3 are not optimal values as the constraint of equality 
does not hold good and also the^pst functions of the plants are non-linear and they are approximated as 
either linear or quadratic. 




It is observed that for anflcSkV GA approach provides optimal values as the constraint of equality is almost 
zero or nearer to zert^Orfinimize the objective function and there is no need of approximating the cost 
functions. Hence ^rtSuWre accurate method and may replace effectively the conventional methods. 

The cost func*for\of two thermal plants of taken as same as in the case of three thermal plants and the cost 
function o^Vnfcined cycle co-generation plant (CCCP) can be taken as F3=(o. 003210) P3 A 2+(7-927738) 
P3+56.6f^^4 Rs/h. The loss co-efficients, operating ranges and the other data used in GA are same as in 
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V. Conclusions 

The genetic algorithm is a searching or optimizing algorithm 



Tables. 5.1 and 5.2 * 

As combined cycle co-generation plant (CCCP) is having low cost, total cost of two thermaJfrJ^jIs with one 
combined cycle co-generation plant is lower than the total cost of three thermal planf^ystVn for any load 
and it is proved by the GA approach given in the table 5.2. The method is accufc^l^md may replays 
effectively the conventional practices presently performed in different central loaidi^j^ched centers. This 
approach is only viable solution for power system with combined cycle co-geo^l^bn plants and it is not 
possible to solve such a ELD problem by conventional technique and it is giverf pVe 5.2 

ba^oV^B the natural evolution principles. 
Because of its capability to solve nonlinear optimization problert^Ww application of GAs to power system 
is a promising area to explore. , 

There is no need of approximating the cost functions oflffh^plants in the differentiable form or quadratic 
form as in the case of conventional methods as Gaj^l* Algorithms are used for solving the functions 
whether they are differentiable or not, continuous*s\not effectively. They do not stuck into local optima, 
because these begin with many initial points andN^arch for the most optimum in parallel and these have 
been proved to be effective and quite robust^r^olving the optimization problems and can handle discrete 
control variable effectively and provide glof^^^lutions. 

VwLscope for Enhancement 

Movement in the simple genetp^lgonthm is accomplished using three primary operators: reproduction, 



crossover and mutation. A gi 



Refined genetic algorith^ 
faster convergence. T^j 




Algorithm works with a population of chromosomes. 



fer from simple genetic algorithms by some improvements made to ensure 
Je elitism, changing probabilities of mutation and crossover.. 



iety of current applications in science, engineering business and the social sciences. 
10ns like multi-objective optimization etc. the advanced operators like Dominance, 
ance, Inversion and other Recording operators may be used. 
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